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Abstract 

In a number of physically important cases, the nonholonomically 
(nonintegrable) constrained Ricci flows can be modelled by exact solu- 
tions of Einstein equations with nonhomogeneous (anisotropic) cosmo- 
logical constants. We develop two geometric methods for constructing 
such solutions: The first approach applies the formalism of nonholo- 
nomic frame deformations when the gravitational evolution and field 
equations transform into systems of nonlinear partial differential equa- 
tions which can be integrated in general form. The second approach 
develops a general scheme when one (two) parameter families of exact 
solutions are defined by any source-free solutions of Einstein's equa- 
tions with one (two) Killing vector field(s). A successive iteration pro- 
cedure results in a class of solutions characterized by an infinite number 
of parameters for a non-Abelian group involving arbitrary functions 
on one variable. We also consider nonlinear superpositions of some 
mentioned classes of solutions in order to construct more general in- 
tegral varieties of the Ricci flow and Einstein equations depending on 
infinite number of parameters and three/ four coordinates on four/ five 
dimensional (semi) Riemannian spaces. 
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1 Introduction 



In recent years, much work has been done on Ricci flow theory and 
fundamental problems in mathematics (USE] (see [H El E] for reviews 
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and references therein). In this context, a number of possible applications 
in modern gravity and mathematical physics were proposed, for instance, 
for low dimensional systems and gravity [HI [9l [10] and black holes and 
cosmology [11] [E]. Such special cases were investigated following certain 
low dimensional or approximative solutions of the evolution equations. 

There were also examined possible connections between Ricci flows, soli- 
tonic configurations and Einstein spaces |13 t ll4 ^ [T5]. In our works, we tackled 
the problem of constructing exact solutions in Ricci flow and gravity theo- 
ries in a new way. Working with general (pseudo) Riemannian spaces and 
moving frames, we applied certain methods from the geometry of Finsler— 
Lagrange spaces and nonholonomic manifolds provided with nonlinear con- 
nection structure (N-connection) [HI [T71 [18l [19] . 

Prescribing on a manifold some preferred systems of reference and sym- 
metries, it is equivalent to define some nonintegrable (nonholonomic, equiv- 
alently, anholonomic) distributions with associated N-connections. From 
this point of view, of the geometry of so-called nonholonomic manifolds, it 
is possible to elaborate a unified formalism for locally fibred manifolds and 
vector /tangent bundles when the geometric constructions are adapted to the 
N-connection structure. We can consider different classes of metric and N- 
connection ansatz and model, for instance, a Finsler, or Lagrange, geometry 
in a (semi) Riemannian (in particular, Einstein) space. Inversely, we can de- 
fine some effective Lagrangian, or Finsler like, fundamental functions for lifts 
of geometric objects for a theory of gravity to tangent bundles in order to 
elaborate a geometric mechanics model for such gravitational and/or gauge 
field interactions, see examples and details in Refs. [HI [201 [ZD [HI [23]. It 
was also proved that constraining some classes of Ricci flows of (semi) Rie- 
mannian metrics we can model Finsler like geometries and, inversely, we 
can transform Finsler-Lagrange metrics and connections into Riemannian, 
or Riemann-Cartan ones [13 ^ 114 ^ 12^. 

The most important idea in constructing exact solutions by geometric 
methods is that we can consider such nonholonomic deformations of the 
frame and connection structures when the Cartan structure equations, Ricci 
flow and/or Einstein equations transform into systems of partial differential 
equations which can be integrated in general form, or one can be derived 
certain bi-Hamilton and solitonic equations with corresponding hierarchies 
and conservation laws, see [191 [E] and references therein. 

The first examples of physically valuable exact solutions of nonholonomic 
Ricci flow evolution equations and gravitational field equations were con- 
structed following the so-called anholonomic frame method [25[ 1261127]. We 
analyzed two general classes of solutions of evolution equations on time like 
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and/or extra dimension coordinate (having certain nontrivial limits to exact 
solutions in gravity theories): The first class was elaborated for solitonic and 
pp-wave nonholonomic configurations. The second class was connected to a 
study of nonholonomic Ricci flow evolutions of three and four dimensional 
(in brief, 3D and 4D) Taub-NUT metrics. Following those constructions 
and further geometric developments in Refs. [13\ [T3j, we concluded that 
a number of important for physical considerations solutions of Ricci flow 
equations can be defined by nonholonomically generalized Einstein spaces 
with effective cosmological constant running on evolution parameter, or (for 
more general and/or normalized evolution flows) by 'nonhomogeneous' (lo- 
cally anisotropic) cosmological constants. 

This is the forth paper in a series of works on nonholonomic Ricci flows 
modelled by nonintegrable constraints on the frame structure and evolu- 
tion of metrics \13\ HM I15j . It is devoted to geometric methods of con- 
structing generic off-diagonal exact solutions in gravity and Ricci flow the- 
ory0 The goal is to elaborate a general scheme when starting with certain 
classes of metrics, frames and connections new types of exact solutions are 
constructed following some methods from nonholonomic spaces geometry 
[201 ED EH [231 IB] and certain group ideas [3ll [32]. The approach to gen- 
erating vacuum Einstein metrics by parametric nonholonomic transforms 
was recently formulated in Ref. [33] (this article proposes a "Ricci flow de- 
velopment" of sections 2 and 3 in that paper). Such results seem to have 
applications in modern gravity and nonlinear physics: In the fifth part- 
ner paper [28], we show how nonholonomic Ricci flow evolution scenaria of 
physically valuable metrics can be modelled by parametric deformations of 
solitonic pp-waves and Schwarzschild solutions. 

One should be noted that even there were found a large number of exact 
solutions in different models of gravity theory [291 [301 [T9l [20l [211 [221 [23] . 
and in certain cases in the Ricci flow theory [71[8l[9l[Tn l [25 l [26 l [27 1 [T5]. 
one has been elaborated only a few general methods for generating new 
physical solutions from a given metric describing a real physical situation. 
For quantum fields, there were formulated some approximated approaches 
when (for instance, by using Feynman diagrams, the formalism of Green's 
functions, or quantum integrals) the solutions are constructed to represent 
a linear or nonlinear prescribed physical situation. Perhaps it is unlikely 
that similar computation techniques can be elaborated in general form in 
gravity theories and related evolution equations. Nevertheless, certain new 

^We shall follow the conventions from the first two partner works in the series; the 
reader is recommended to study them in advance. 
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possibilities seem to be opened after formulation of the anholonomic frame 
method with parametric deformations for the Ricci flow theory. Although 
many of the solutions resulting from such methods have no obvious physical 
interpretation, one can be formulated some criteria selecting explicit classes 
of solutions with prescribed symmetries and physical properties. 

The paper has the following structure: In section 2, we outline some re- 
sults on nonholonomic manifolds and Ricci flows. Section 3 is devoted to the 
anholonomic frame method for constructing exact solutions of Einstein and 
Ricci flow equations. There are analyzed the conditions when such solutions 
define four and five dimensional foliations related to Einstein spaces and 
Ricci flows for the canonical distinguished connection and the Levi Civita 
connection. In section 4, we consider how various classes of metrics can be 
subjected to nonholonomic deformations and multi-parametric transforms 
(with Killing symmetries) resulting in new classes of solutions of the Ein- 
stein/ Ricci flow equations. We consider different ansatz for metrics and 
two examples with multi-parametric families of Einstein spaces and related 
Ricci flow evolution models. The results are discussed in section 5. The 
reader is suggested to see Appendices before starting the main part of the 
paper: Appendix A outlines the geometry of nonlinear connections and the 
anholonomic frame method of constructing exact solutions. Appendix B 
summarizes some results on the parametric (Geroch) transforms of vacuum 
Einstein equations. 

Notation remarks: It is convenient to use in parallel two types of de- 
notations for the geometric objects subjected to Ricci flows by introducing 
"left-up" labels like -^7 = 7(...,x). Different left-up labels will be also con- 
sidered for some classes of metrics defining Einstein spaces, vacuum solutions 
and so on. We shall also write "boldface" symbols for geometric objects and 
spaces adapted to a noholonomic (N-connection) structure, for instance, 
V, E, ... A nonholonomic distribution with associated N-connection struc- 
ture splits the manifolds into conventional horizontal (h) and vertical (v) 
subspaccs. The geometric objects, for instance, a vector X can be written 
in abstract form X = (hX,vX), or in coefficient forms, X" = = 
X-) , equivalently decomposed with respect to a general nonholonomic 
frame = (ej,ea) or coordinate frame da = (dyda) for local h- and v- 
coordinates u = {x,y), or u" = when da = d/du— and di = d/dx-, 

da = d/dx-. The h-indices i, j, fc, ... = 1, 2, ...n will be used for nonholonomic 
vector components and the v-indices a,b,c... = n + l,n + 2, ...n + m will be 
used for holonomic vector components. Greek indices of type a,/?,... will 
be used as cumulative ones. We shall omit labels, indices and parametric/ 
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coordinate dependencies if it does not result in ambiguities. 



2 Preliminaries 

In this section we present some results on nonholonomic manifolds and 
Ricci flows [131 [T3] selected with the aim to outline a new geometric method 
of constructing exact solutions. The anholonomic frame method and the 
geometry of nonlinear connections (N-connections) are considered, in brief, 
in AppendixlAl The ideas on generating new solutions from one/ two Killing 
vacuum Einstein spacetimes [311 [32] are summarized in Appendix |Bl 

2.1 Nonholonomic (pseudo) Riemannian spaces 

We consider a spacetime as a (necessary smooth class) manifold V of 
dimension n + m, when n > 2 and m > 1 (a splitting of dimensions being 
defined by a N-connection structure, see ()A.3P ). Such manifolds (equiva- 
lently, spaces) are provided with a metric, g = gape°' e", of any (pseudo) 
Euclidean signature and a linear connection D = {T^^^e'^} satisfying the 
metric compatibility condition Dg = oH The components of geometrical 
objects, for instance, gap and T"^^, are defined with respect to a local base 

(frame) and its dual base (co-base, or co-frame) e" for which e^J = , 
where "J" denotes the interior product induced by g and 5a is the Kronecker 
symbol. For a local system of coordinates ti" = (x*,y") on V (in brief, 
u = {x,y)), we write respectively 

d d 

ea = (ci = di = ea = da = ^) and e'^ = {e^ = dx^ = dy^), 

for eaje"^ = 5^; the indices run correspondingly values of type: = 
1, 2, n and a, b, ... = n + 1, n + 2, n + m for any conventional splitting 
a = {i,a),(3 = {j,b), ... 

Any local (vector) basis and dual basis can be decomposed with 
respect to other local bases Ca' and 6^ by considering frame transforms, 

ea = A^' {u)ea' and = (n)e^' , (1) 

where the matrix A^^, is the inverse to Aa' . It should be noted that an 
arbitrary basis is nonholonomic (equivalently, anholonomic) because, in 

■^in this work, the Einstein's summation rule on repeating "upper-lower" indices will 
be applied if the contrary will not be stated 
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general, it satisfies certain anholonomy conditions 

eaef3 - epea = Wj^g (2) 

with nontrivial anholonomy coefficients ^^^(u). For VFj^ = 0, we get holo- 
nomic frames: for instance, if we fix a local coordinate basis, = da- 

Denoting the covariant derivative along a vector field X = X'^Ca as 
Dx = X\D, we can define the torsion 

T(X, Y) = DxY- DyX - [X, Y] , (3) 

and the curvature 

7^(X, Y)Z = DxDyZ - DyDxZ - D[x,y]Z, (4) 

tensors of connection D, where we use "by definition" symbol "=" and 
[X,Y] = XY - YX. The components T = {T°^} and TZ = are 
computed by introducing X— >eQ,y^e/3,Z— into respective formulas 
and (H). 

The Ricci tensor is constructed Ric{D) = {Rp^y = R'^p^a}- "^^^ scalar 
curvature R is by definition the contraction with g"^^ (being the inverse 
to the matrix gap), R =1= g'^^Raf3-, and the Einstein tensor is £" = {Eap == 
Rajd — ^dajsR}- The vacuum (source-free) Einstein equations are 

£ = {Eap = Ra^} = 0. (5) 

In general relativity theory, one chooses a connection D = V which is 
uniquely defined by the coefficients of a metric, gap, following the conditions 
of metric compatibility, Vg' = 0, and of zero torsion, ,T = 0. This is the 
so-called Levi Civita connection ,D = V. We shall respectively label its 
curvature tensor, Ricci tensor, scalar curvature and Einstein tensor in the 
form ,7^ = { ,Ric{V) = { ,Rap = .R^'p^a}^ ^ i^"/3 and 

i-^^ = { lEap}. 

Modern gravity theories consider extra dimensions and connections with 
nontrivial torsion. For instance, in string gravity [35\ I36j. the torsion coef- 
ficients are induced by the so-called anti-symmetric iif-fields and contain 
additional information about additional interactions in low-energy string 
limit. A more special class of gravity interactions are those with effective 
torsion when such fields are induced as a nonholonomic frame effect in a 
unique form ^ by prescribing a nonholonomic distribution (|A.3p . defining 
a nonlinear connection structure, N-connection, on a (pseudo) Riemannian 
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manifold V enabled with a metric structure ()A.9p . Such spaces with lo- 
cal fibred structure are called nonholonomic (in more special cases, when 
the nonholonomy is defined by a N-connection structures, the manifolds are 
called N-anholonomic) |171ll9j. The N-anholonomic spaces can be described 
in equivalent form by two linear connections V ()A.16p and D ()A.17p . both 
metric compatible and completely stated by a metric ()A.10p . equivalently 
As a matter of principle, the general relativity theory can be for- 
mulated in terms of both connections, V and D; the last variant being with 
nonholonomic constraints on geometrical objects. One must be emphasized 
that the standard approach follows the formulation of gravitational field 
equations just for the Einstein tensor ^£ for V which, in general, is different 
from the Einstein tensor E for D. 

A surprising thing found in our works is that for certain classes of generic 
off— diagonal metric ansatz ()A.9P it is possible to construct exact solutions in 
general form by using the connection D but not the connection V. Here we 
note that having defined certain integral varieties for a first class of linear 
connections we can impose some additional constraints and generate solu- 
tions for a class of Levi Civita connections, for instance, in the Einstein 
and string, or Finsler like, generalizations of gravity. Following a geometric 
N-adapted formalism (the so-called anholonomic frame method), such solu- 
tions were constructed and studied in effective noncommutative gravity [18] , 
various locally anisotropic (Finsler like and more general ones) extensions of 
the Einstein and Kaluza-Klein theory, in string an brane gravity [21} [23] 
and for Lagrange-Fedosov manifolds [34], see a summary in [19]. 

The anholonomic frame method also allows us to construct exact solu- 
tions in general relativity: One defines a more general class of solutions for 
D and then imposes certain subclasses of nonholonomic constraints when 
such solutions solve the four dimensional Einstein equations for V. Here we 
note that by nonholonomic deformations we were able to study nonholo- 
nomic Ricci flows of certain classes of physically valuable exact solutions 
like solitonic pp-waves [25] and Taub NUT spaces [261 127j . In this work, we 
develop the approach by applying new group methods. 

Certain nontrivial limits to the vacuum Einstein gravity can be selected 
if we impose on the nonholonomic structure such constraints when 

^= (6) 

even, in general, D 7^ V. We shall consider such conditions when D and V 

^In this work, we shall use only V and the canonical d-connection ID and, for simplicity, 
we shall omit "hat" writing D if that will not result in ambiguities. 
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have the same components with respect to certain preferred bases and the 
equahty can be satisfied for some very general classes of metric ansatz0 
We shall use left-up labels "o" or for a metric, 

= °gap ® or ^5 = ^g^p e" 

being (correspondingly) a solution of the vacuum Einstein (or with cosmo- 
logical constant) equations £^ = ([5]) or of the Einstein equations with a 
cosmological constant A, R^p = \gap-, for a linear connection D with possi- 
ble torsion T 7^ 0. In order to emphasize that a metric is a solution of the 
vacuum Einstein equations, in any dimension n + m > 3, for the Levi Civita 
connection V, we shall write 

°g = °gaf3 ef" or = ^g^p ® e\ 

where the left-low label will distinguish the geometric objects for the 
Ricci flat space defined by a Levi Civita connection V. 

Finally, in this section, we note that we shall use "boldface" symbols, 
for instance, if = ^g^/j e" ® defines a nonholonomic Einstein space 
as a solution of 

R-a/? = Aga/3 (7) 

for the canonical d-connection D. 

2.2 Evolution equations for nonholonomic Ricci flows 

The normalized (holonomic) Ricci flows O [H O [6] for a family of metrics 
go/six) = gaf3iu-,x)j parametrized by a real parameter with respect to 
the coordinate base da = d/du—, are described by the equations 

d 2r 

^^o^ = -2 ,Raj3 + -^ga/3, (8) 

where the normalizing factor r = J ,RdV/dV is introduced in order to 
preserve the volume V. H 



*We emphasize that difFerent hnear connections may be subjected to difTerent rules of 
frame and coordinate transforms. It should be noted here that tensors and nonlinear and 
linear connections transform in different ways under frame and coordinate changing on 
manifolds with locally fibred structures, see detailed discussions in |16l 119) . 

^The Ricci flow evolution equations were introduced by R. Hamilton [l], as evolution 
equations 

99 M 



--^ = -2 &,{X), 
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For N-anholonomic Ricci flows, the coefficients ggjs are parametrized in 
the form (|A.9p . see proofs and discussion in Refs. [131 113 113 HZ!- With 
respect to N-adapted frames ()A.4p and (jA.Sp . the Ricci flow equations ([8|), 
redefined for V ^ D and, respectively, ,Raf3 R-a/? are 



d_ 
d_ 

t; 9ab 
ox 



and Rai = 0, 



Rii + >^9; 



d_ 

'dx^ 



g,a^{N^Nf), (9) 

(10) 
(11) 



where A = r/5 the Ricci coefficients R^^ and are computed with respect 
to coordinate coframes. The equations ([lip constrain the nonholonomic 
Ricci flows to result in symmetric metricsO 

Nonholonomic deformations of geometric objects (and related systems of 
equations) on a N-anholonomic manifold V are defined for the same metric 
structure g by a set of transforms of arbitrary frames into N-adapted ones 
and of the Levi Civita connection V into the canonical d-connection D, 
locally parametrized in the form 



da = idi,da) ^ea = (ei,ea); gaj3 [gij, gab, N"-]; ,r' 



al3 



A rigorous proof for nonholonomic evolution equations is possible following 
a N-adapted variational calculus for the Perelman's functionals presented 
in Refs. [U]. For a five dimensional space with diagonal d-metric ansatz 
(jA.lOp . when gij = diag[iil,g2-,gz\ and gab = diag[g4,, g^], we considered [25j 
the nonholonomic evolution equations 



d_ 
d} 
d_ 

dx 
R 



o 9i: 

OX 



9aa 



Rii '^9ii 



~2 yRaa ~ X9aa j j 

for a ^ f3, 



(12) 

(13) 
(14) 

parametrized by a real x shall underline symbols or indices in order to emphasize 
that certain geometric objects/ equations are given with the components defined with 
respect to a coordinate basis). For our further purposes, on generalized Riemann-Finsler 
spaces, it is convenient to use a different system of denotations than those considered by 
R. Hamilton or Grisha Perelman on holonomic Riemannian spaces. 

®In Refs. [131 124) . we discuss this problem related to the fact that the tensor is 
not symmetric which results, in general, in Ricci flows of nonsymmetric metrics. 



for a set of Riemannian metrics g ^(x) ^^i^d corresponding Ricci tensors 
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with the coefficients defined with respect to N-adapted frames (|A.4p and 
(jA.Sp . This system can be transformed into a similar one, like ([9|)- pip . by 
nonholonomic deformations. 



3 Off- Diagonal Exact Solutions 

We consider a five dimensional (5D) manifold V of necessary smooth 
class and conventional splitting of dimensions dimV = n + m for n = 3 
and m = 2. The local coordinates are labelled in the form = {x'^,y"') = 
(x^, X*, = V, y^), for i = 1, 2, 3 and i = 2, 3 and a, b, ... = 4, 5. Any coordi- 
nates from a set can be a three dimensional (3D) space, time, or extra 
dimension (5th) one. Ricci flows of geometric objects will be parametrized 
by a real x- 

3.1 Off— diagonal ansatz for Einstein spaces and Ricci flows 

The ansatz of type (jA.lOp is parametrized in the form 

+h4(^x^,v^ 5v ^ 6v + (^x'' , 5y<^5y, 
5v = dv + Wi (^x^,v^ dx\ 6y = dy + rii dx^ (15) 

with the coefficients defined by some necessary smooth class functions 

gi = ±1,52,3 = 52,3(2;^, a;^),/i4,5 = /i4,5 (a;*,?;), 
Wi = Wi{x\v),ni = ni{x\v). 

The off-diagonal terms of this metric, written with respect to the coordi- 
nate dual frame tin" = (dx*,(iy°), can be redefined to state a N-connection 
structure N = [N^ = Wi{x'' ,v),N-' = ni{x^,v)] with a N-elongated co-frame 
()A.5p parametrized as 

= dx'^, (? = dx"^, e' = dx^, 
e'^ = 6v = dv + Widx^, = 6y = dy + riidx^ . (16) 

This coframe is dual to the local basis 

d f k \ 9 f k \ d d d 

ei = — -Wi[x ,v] ni[x ,v] -^,64 = -^,65 = (17) 

ax* V J ov \ J oy° ov oy° 
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We emphasize that the metric (jlSp does not depend on variable , i.e. it 
posses a KiUing vector 65 = d/dy^ ^ and distinguishes the dependence on the 
so-cahed "anisotropic" variable y'^ = v. 

The above considered ansatz and formulas can be generalized in order 
to model Ricci flows, 

g = gidx dx + g2{x ,x ,x)dx ® dx + 53 (x ,x ,x)dx <^dx 

^5v = dv + Wi (^x'' ,v,xj dx\ ^5y = dy + Ui ^x^ ,v,xj dx^ (18) 
with corresponding flows for N-adapted bases, 

ea = (ei,ea) ^ = ( ^ei,ea) = e„(x) = (ei(x),ea), 
e° = (e\e'^) ^ = {e\ ^e'') = e"(x) = (e*,e"(x)) 

defined by Wi [x^,v) — > Wi {x^ ,v, A) , rii {x^ ,v) Ui (x'^, v, A) in (dZl), 

Computing the components of the Ricci and Einstein tensors for the 
metric (|18p (see main formulas in Appendix and details on tensors compo- 
nents' calculus in Refs. [18\ I19j). one proves that the corresponding family 
of Rice tensors for the canonical d-connection with respect to N-adapted 
frames are compatible with the sources (they can be any matter fields, string 
corrections, Ricci flow parameter derivatives of metric, ...) 

= [Tl = T2 + T4,Ti = T2(x2,x3,t;,x),Ti = T2(x2,x3,t-,x), 

= T4(x^x^x),Ti = T4(x^x^x)] (19) 

transform into this system of partial differential equations: 
R| = l4(x) (20) 

_ 1 .*-95_^M)!_,.._^M3_^(^_^;',,_T^,,._,,_^)_ 



25253 2^2 253 253 2^2 

St = Slix) = [K (invlW)* -/^r] = -T2(x2,x3,^,x), (21) 

P / X /%) Q»(X) n 

= -"'•<'='2MJ)-2MJ)=°' '''' 

ftsi = -|^["r(O+7(xK(Ol = 0. (23) 
where, for /i| 5 / 0, 

ai(x-) = hHx)a,4,(x). m = K(X) 4>'(X). (24) 

2/i5(x) /i4(x) ^Mx)Mx)l 
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when the necessary partial derivatives are written in the form a' = da/dx^ , 
a' = da/dx^, a* = da/dv. In the vacuum case, we must consider T2,4 = 0. 
We note that we use a source of type (|19p in order to show that the anholo- 
nomic frame method can be apphed also for non-vacuum configurations, 
for instance, when T2 = A2 = const and T4 = A4 = const, defining local 
anisotropics generated by an anisotropic cosmological constant, which in its 
turn, can be induced by certain ansatz for the so-called H-iield (absolutely 
antisymmetric third rank tensor) in string theory [18\ I19j . We note that 
the off-diagonal gravitational interactions and Ricci flows can model locally 
anisotropic configurations even if A2 = A4, or both values vanish. 

Summarizing the results for an ansatz (jlSp with arbitrary signatures 
Co = (ei, £2, £3, £4, £5) (where ea = ±1) and hi ^ and /15 7^ 0, for a fixed 
value of X, one proves, see details in [I8l[l9], that any off — diagonal metric 



°g = ei dx'^ (gi dx^ + €2g2ix'') dx'^ (gi dx'^ + e3g3{x'') dx^ dx^ 

+e4hl{x') [f* {x\ v)]^ I? {x\ v)\5v(^ 5v 

+65 [f {x\v) - foix' f 6y'>®6y\ 
Sv = dv + Wk {x\v) dx^, 6y^ = dy^ + n^ {x\ v) dx^ , (26) 

with the coefficients being of necessary smooth class and the idices with 
"hat" running the values = 2,3, where g-j^ ^x'^j is a solution of the 

2D equation (pOj) for a given source T4 ( x* 



, {x\v) = ?[oi (xO - jhlix^) I T2{x^,v)r {x\v) [f {x\v) - /o(x*)] dv, 

and the N-connection coefficients Nf = Wi{x^,v), Nf = ni{x^,v) are com- 
puted following the formulas 

nk = (x') + nfc[2] (x^) j ^[ ^•^'''^) <;{x\v)dv, (28) 

J [/(x%w) - /o(x*)] 

define an exact solution of the system of Einstein equations d?]) . It should be 
emphasized that such solutions depend on arbitrary functions / (x*, v) , with 

/* / 0, /o(x'), /lo(x'), ?[o] {x') (x*) , nfc[2] (x*) and T2(x^, w), T4 (xM . 
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Such values for the corresponding signatures = ±1 have to be stated by 
certain boundary conditions fohowing some physical considerationslll 

The ansatz of type ([15]) with /14 = but 7^ (or, inversely, ^ but 
= 0) consist more special cases and request a bit different methods for 
constructing exact solutions. Nevertheless, such solutions are also generic 
off— diagonal and they may be of substantial interest (the length of paper 
does not allow us to include an analysis of such particular cases) . 

3.2 Generalization of solutions for Ricci flows 

For families of solutions parametrized by Xi consider flows of the 
generating functions, g2{x'^,x)j or 53(x*,x)) and /(x*,f,x), and various 
types of integration functions and sources, for instance, n;j[i] (x*,x) and 

n^,[2] {x^,x) and T2{x^ ,v,x), respectively, in formulas ([27j) and (p8]) . Let us 
analyze an example of exact solutions of equations ()12p - (ll4p : 
We search a class of solutions of with 

92 = e2Tu{x'^,x^,x),93 = e3w{x'^,x^,x), 

for a family of ansatz (llSp with any prescribed signatures = ±1 and 
non-negative functions zu and h. Following a tensor calculus, adapted to 
the N-connection, for the canonical d-connectionH we express the integral 
variety for a class of nonholonomic Ricci flows as 

e2(ln|w|)" + e3(ln|zu|)" = IX-h^d^iuif, (29) 
hi = /i<;4 

for 

?4(x ,x ,v) = ?4[o](x ,x ) - J J -f^dv 

^\ = /t[o](a;') (Vl^5 {x\x\v)\y (30) 

''Our classes of solutions depending on integration functions are more general than 
those for diagonal ansatz depending, for instance, on one radial like variable like in the 
case of the Schwarzschild solution (when the Einstein equations are reduced to an effective 
nonlinear ordinary differential equation, ODE). In the case of ODE, the integral varieties 
depend on integration constants to be defined from certain boundary/ asymptotic and 
symmetry conditions, for instance, from the constraint that far away from the horizon the 
Schwarzschild metric contains corrections from the Newton potential. Because our ansatz 
(|15|) transforms ((TJ in a system of nonlinear partial differential equations transforms, the 
solutions depend not on integration constants but on integration functions. 

**similar computations are given in [TS] and Chapter 10 of [19) 
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and, for ip = — In 



VWJ^\/\K\ 



m = 0,n2 = ns = n[i]{x'^,x^,x) +n[2]ix'^,x^,x) j dv h4^/ (^^/]h^^ , 

where the partial derivatives are denoted in the form ip' = dip/dx'^,ip' = 
dif/dx'^, if* = dif/dv, = d/dx'^, and arbitrary /15 when ^ 0. For A = 0, 
we shah consider <;"4[o] = 1 and /i[o](x*) = const in order to solve the vacuum 
Einstein equations. There is a class of solutions when 



for a function C(x^,x^). This is compatible with the condition (|30p . and 
we can chose such configurations, for instance, with njij = and any 
n[2]{x'^ , x^ , x) and ■uj{x'^,x^,x) solving the equation ([29]) . 

Putting together (f29]) - (f3T] ). we get a class of solutions of the system (fT2]) - 
(fHI) (the equations being expressed equivalently in the form (p0l) - ([23]l ) for 
nonholonomomic Ricci flows of metrics of type (jlSp . 

'^g = eidx"^ (8) dx"^ + tz7(x^, x^, x) [e2(ix^ (g) dx^ + esdx^ (g) dx^] 

+/i4 (x^, x^, 6v 6v + [x'^ ,x^ ,v) ^5y(^ ^6y, 
6v = dv + W2{x'^ iX^ ,v) dx^ + w^{^x^ ,x^ ^v) dx^ ^ (32) 
^6y = dy + n2 [x"^ iX^ ,v,x) [dx^ + dx^]. 

Such solutions describe in general form the Ricci flows of nonholonomic 
Einstein spaces constrained to relate in a mutually compatible form the 
evolution of horizontal part of metric, w{x'^ ^x^ ,x)i with the evolution of 
N-connection coefficients f^2 = "-3 = n2 (x^,x^,f,x) . We have to impose 
certain boundary/ initial conditions for x = 0, beginning with an explicit 
solution of the Einstein equations, in order to define the integration functions 
and state an evolution scenario for such classes of metrics and connections. 



3.3 4D and 5D Einstein foliations and Ricci flows 

The method of constructing 5D solutions can be restricted to generate 
4D nonholonomic configurations and generic off-diagonal solutions in gen- 
eral relativity. In order to consider reductions bD AD for the ansatz (jl5p . 
we can eliminate from the formulas the variable x^ and consider a 4D space 

(parametrized by local coordinates (x^, x^, w, y^)) trivially embedded 
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into a 5D spacetime V (parametrized by local coordinates (^x^ , x'^ , , v , y^) 
with gii = ±1,(713 = 0,a = 2,3,4,5). In this case, there are possible 
4D conformal and anholonomic transforms depending only on variables 
of a 4D metric g^j^ (^x'^ , x^ , of arbitrary signature. To em- 
phasize that some coordinates are stated just for a 4D space we might use 
"hats" on the Greek indices, 3, /?, ... and on the Latin indices from the mid- 
dle of the alphabet, i, j, ... = 2,3; local coordinates on are parametrized 
u°' = ^x*, y'^^ = (x^, x^, = V, y^) , for a, b, ... = 4, 5. The ansatz 

g = 52 dx'^ (g) dx^ + gs dx^ ® dx^ + 5v (g) 6v + 5y^ 5y^, (33) 
is written with respect to the anholonomic co-frame (^dx^ , Sv , 6y^^ , where 

6v = dv + vrjdx^ and 5y^ = dy^ + ■nqdx'' (34) 
is the dual of (5^, d^, d^) , for 

&>=&■+ ur.d4 + n^5, (35) 
and the coefficients are necessary smoothly class functions of type: 

9j = 9j{x''),h4^5 = h4^5{x'',v), 

nr. = W2{x^,v),-nq = 'nq{x^,v)] i,k = 2,2>. 

In the 4D case, a source of type ()19p should be considered without the 
component T} in the form 

T| = diag[Tl = Ti = T2(x^^;), T| = = T4(x^)]. (36) 

The Einstein equations with sources of type (j36p for the canonical d-connec- 
tion ()A.16P defined by the ansatz (j33p transform into a system of nonlinear 
partial differential equations very similar to (j20p - (j23p . The difference for 
the 4D equations is that the coordinate x^ is not contained into the equa- 
tions and that the indices of type .. = 1,2,3 must be changed into the 
corresponding indices z, j, .. = 2, 3. The generated classes of 4D solutions are 
defined almost by the same formulas (p6|) . ([271) ^"^^ (pHj) . 

Now we describe how the coefficients of an ansatz (|33p defining an exact 
vacuum solution for a canonical d-connection can be constrained to generate 
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a vacuum solution in Einstein gravity: We start with the conditions (|A.20p 
written (for our ansatz) in the form 

^-w^hl-2wlhi = 0, (37) 

= 0, (38) 

ox" 

n^h^ = 0. (39) 

These equations for nontrivial values of w-^ and constructed for some 
defined values of and /15 must be compatible with the equations (|2ip - 
()23p for T2 = 0. One can be taken nonzero values for w-^ in \22\\ if and only 
if = because the the equation ([2T|) imposes the condition /3 = 0. This 
is possible, for the sourceless case and /15 7^ 0, if and only if 



(/) = In h\/ ■\/\hih^\ = const ^ (40) 

see formula ([25]) . A very general class of solutions of equations ([37|) . ([38|) 
and (I40p can be represented in the form 

/i4 = e^hl (6*)^/i5 = €5(6 + 60)', (41) 
9(6 + 60) 



where ho = const and 6 = b{x^,v) is any function for which 6* ^ and 
60 = 6o(x'^) is an arbitrary integration function. 

The next step is to satisfy the integrability conditions (lA.lSp defin- 
ing a foliated spacetimes provided with metric and N-connection and d- 
connection structures [18^ \T9\ [34] (we note that (pseudo) Riemannian foli- 
ations are considered in a different manner in Ref. [ITJ) for the so-called 
Schouten - Van Kampen and Vranceanu connections not subjected to the 
condition to generate Einstein spaces). It is very easy to show that there 
are nontrivial solutions of the constraints (|A.18P which for the ansatz (j33p 
are written in the form 

■^^2 - ^3 + ^3'«^2 - = 0, (42) 

^2 ~ + ^3^2 ~ ^2f^3 = 0. 

We solve these equations for rig = ^3 = if we take any two functions 
"-2,3(2;^) satisfying 

n'2 -nl = (43) 
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(this is possible by a particular class of integration functions in (|28p when 

n^j2] {^^^ ~ ^ '^fcfi] (■^*) constraint to satisfy just the conditions 

(j43|) l. Then we can consider any function b{x^,v) for which w-j^ = (b*)^^ 
d^{b + 60) solve the equation (jl2]) . In a more particular case, one can be 
constructed solutions for any v), b* 7^ 0, and n2 = and na = ns{x^,v) 
(or, inversely, for any n2 = n2{x'^,v) and 723 = 0). We also note that the 
conditions (|A.19p are solved in a straightforward form by the ansatz p3p . 
We conclude that for any sets of coefficients 

h4{x'',v), h^{x^,v), ^^;r(x^^;), n2,z{x'") 



respectively generated by functions b{x ,v) and n^j-^j \^x^j , see (j^T]) . and 

satisfying (jl3]) . the generic off-diagonal metric ([33|) possess the same coef- 
ficients both for the Levi Civita and canonical d-connection being satisfied 
the conditions ([6]) of equality of the Einstein tensors. Here we note that any 
2D metric can be written in a conformally flat form, i. e. we can chose such 
local coordinates when 



92{dx^f + gsidx^f = e^(^') e^{dx'Y + e^{dx 



„2^2 , , fj^3\2 



for signatures = ±1, in (f33]) . 

Summarizing the results of this section, we can write down the generic 
off-diagonal metric (it is a 4D dimensional reduction of (j26p ) 



og ^ ,^') dx^ ® dx^ + eg (ix^ dx^] (44) 

+64/10 {x\v)Y (g) 6v 
+£5 [b (x^, x^ v) - 6o(x^ x^)] ^ 6y^ (g) 

Ji; = dv + W2 {x'^ , x"^ , v) dx"^ + W3 [x"^ , x"^ , v) dx^ , 
6y^ = dy^ + 77-2 (x^,x^) dx^ + 723 (x^,x^) dx^, 

defining vacuum exact solutions in general relativity if the coefficients are 
restricted to solve the equations 

e2^p" + e3^" = 0, (45) 

w'2 — wl + 'W3W2 — W2VJ^ = 0, 

n'2 - n' = 0, 

for W2 = (b*)"^ (b + 60)' and 7i;3 = (6*)^^ {b + 60)', where, for instance, 
no = 92f^3 and n'o = 93122. 
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We can generalize ()44p similarly to (j26p in order to generate solutions 
for nontrivial sources (j36p . In general, they will contain nontrivial anliolo- 
nomically induced torsions. Such configurations may be restricted to the 
case of Levi Civita connection by solving the constraints (j37p - (|39p in order 
to be compatible with the equations ([21]) and ([^2]) for the coefficients 

and j3 computed for /i^ / and In ht^/ \/\hihh\ = (f>{x'^ , , v) ^ const, 

see formula (|25p . resulting in more general conditions than (j40p and (j4ip . 
Roughly speaking, all such coefficients are generated by any /14 (or /15) de- 
fined from (j22p for prescribed values /15 (or /15) and T2{x^, v). The existence 
of a nontrivial matter source of type (j36p does not change the condition 
ni = 0, see ([39]), necessary for extracting torsionless configurations. This 
mean that we have to consider only trivial solutions of (j23p when two func- 
tions = n^(x^,x^) are subjected to the condition ([42]) . We conclude that 
this class of exact solutions of the Einstein equations with nontrivial sources 
([36]) . in general relativity, is defined by the ansatz 

/i4 (x^,^^,^) (J-u (g) + /i5 (x^,a;^,t>) 6y^(S>6y^, 
5v = dv + W2{x'^ ,x'^ ,v) dx^ + w:i{x^ ^x^ ,v) dx'^ , 
5y^ = dy^ + 722 (x^, x^) dx^ + 723 (x^, x^) dx^, 

where the coefficients satisfy the conditions 

€2^"" + esV'" = ^2 

KHh^h = T2, (47) 

"^2 - ^3 + ^3^2 - ^2-^3 = 0, 

722 -n-l = 0, 

for nr. = d^(j)/(j)*, see ([25]) . being compatible with ([37]) and ([38]) . for given 

sources T4(x'^) and T2(x'^,t;). We emphasize that the second equation in 
([TT]) relates two functions /i4 and . In references [201 ED 1221 1231 HE] , we 
investigated a number of configurations with nontrivial two and three di- 
mensional solitons, reductions to the Riccati or Abbel equation, defining 
off-diagonal deformations of the black hole, wormhole or Taub NUT space- 
times. Those solutions where constructed to be with trivial or nontrivial 
torsions but if the coefficients of the ansatz ([16]) are restricted to satisfy 
the conditions ([47]) in a compatible form with ([37]) and ([38]) . for sure, such 
metrics will solve the Einstein equations for the Levi Civita connection. We 
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emphasize that the ansatz (j46p defines Einstein spaces with a cosmological 
constant A if we put T2 = T4 = A in ([T7|). 

Let us formulate the conditions when famihes of metrics (j46p subjected 
to the conditions (j47p will define exact solutions of the Ricci flows of usual 
Einstein spaces (for the Levi Civita connection). We consider the ansatz 



^g(;^) = e^(^ '^^ [ea dx^ dx'^ + eg dx^ ® dx^] + (48) 

6v = dv + W2 {x'^ , x"^ , v) dx'^ + W3 (^x'^ , x^ , v") dx"^ , 
^5y^ = dy^ + n2 (x^, x^, x) [dx"^ + dx^], 

which is a subfamily of (j32p . when = e^*^^^'^^'-^^ and n2 = na does not 
depend on variable v and the coefficients satisfy the conditions (j29p and 
(j30]) . when np] = but n^ij can be nontrivial in ([3T]) . and (additionally) 

e2V'"(x) + e3/(x) = A, 

hlHh^hr, = A, (49) 

'^2(x) - "-'(x) = 0, 

for nr- = c^(/>/(/>*, see (f25]) . being compatible with (f37|) and (p8]) . for given 
sources T4 = A and T2 = A. The family of metrics (|48p define a self- 
consistent evolution as a class of general solutions of the Ricci flow equations 
(|12p - (|14p transformed equivalently in the form (|2U p - (j23p . The additional 
constraints (09]) define an integral subvariety (foliation) of (j32p when the 
evolution is selected for the Levi Civita connection. 



4 Nonholonomic and Parametric Transforms 

Anholonomic deformations can be defined for any primary metric and 
frame structures on a spacetime V (as a matter of principle, the primary 
metric can be not a solution of the gravitational field equations). Such 
deformations always result in a target spacetime possessing one Killing vec- 
tor symmetry if the last one is constrained to satisfy the vacuum Einstein 
equations for the canonical d-connection, or for the Levi Civita connection. 
For such target spacetimes, we can always apply a parametric transform 
and generate a set of generic off-diagonal solutions labelled by a parame- 
ter 6 (|B.2p . There are possible constructions when the anholonomic frame 
transforms are applied to a family of metrics generated by the parametric 
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method as new exact solutions of the vacuum Einstein equations, but such 
primary metrics have to be parametrized by certain type ansatz admitting 
anholonomic transforms to other classes of exact solutions. Additional con- 
straints and parametrizations are necessary for generating exact solutions of 
holonomic or nonholonomic Ricci flow equations. 

4.1 Deformations and frame parametrizations 

Let us consider a (n + m)-dimensional manifold (spacetime) V, n > 
2,m > 1, enabled with a metric structure g = g (Bn h distinguished in the 
form 

g = giiu)idx'f + haiu){ef, (50) 

The local coordinates are parametrized u = {x,y) = {u" = (x*,y")}, for 
the indices of type i,j,k,... = 1,2, ...,n (in brief, horizontal, or h-indices/ 
components) and a,b,c, ... = n + l,n + 2, ...n + m (vertical, or v-indices/ 
components). We suppose that, in general, the metric ([50]) is not a solution 
of the Einstein equations but can be nonholonomically deformed in order 
to generate exact solutions. The coefficients -/V?(Ti) from ([50j) state a con- 
ventional (n -|- m)-splitting 0^ in any point u € V and define a class of 
'N-adapted' local bases 

and local dual bases (co-frames) c = (c, c), when 

c" = (c> = dx\c^ = dy^ + N^iu) dx^ , (52) 

for cj e = I, i.e. e^J = 5^, where the inner product is denoted by 'J ' and 
the Kronecker symbol is written 5a. The frames (I5ip satisfy the nonholonomy 
(equivalently, anholonomy) relations 

e^e^ - 6/360, = w^^e,^ 

with nontrivial anholonomy coefficients 

^% = -wr,="?,=e,(iVf)-ei(iV;), (53) 
wL = -^l = ea{Nl). 
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A metric g = g (Bn h parametrized in the form 

g = g,{u){df + ga{u){c'^), 
c'^ = dy" + N^{u)dx' 



(54) 



is a nonhlonomic transform (deformation), preserving the (n + m)-sphtting, 
of the metric, g = g h if the coefficients of (j50|) and ([5^ are related by 
formulas 

9i = Vi{u) m, ha = Va{u) ha and iVf = v?{u)K, (55) 

where the summation rule is not considered for the indices of gravitational 
'polarizations' T]a = iVi^Va) and r]f in ([^^ . For nontrivial values of r]f{u), 
the nonholonomic frames ()5ip and (|52p transform correspondingly into 



(56) 



and 



c° = (cJ' = dx*, c" = dy" + TVf (n) dx^) (57) 

with the anholonomy coefficients W^^ defined by ()A.7p . 

We emphasize that in order to generate exact solutions, the gravitational 
'polarizations' r]a = {jli-^Va) and r/f in ([55]) are not arbitrary functions but 
restricted in a such form that the values 



± 1 

2 „3\ 



g2{x ,x 
h/i{x\v) 
Wi{x\v) 



r?4(^x") /i4(n"), h5{x\v) 



(58) 



= r?3(n")53(«"), 

r?5(7X-) /i5(^^"), 

7?f(^.")7Vf(tx"), 



define an ansatz of type (j26p . or ()44p (for vacuum configurations) and (06 
for nontrivial matter sources T2(x^,x^,f) and T4(x^,x^). 
Any nonholonomic deformation 



g = 5®iv ^ — ' g = g®Nh 
can be described by two frame matrices of type (|A.ip . 



5a- 



(59) 



(60) 



generating the d-metric gQ,/3 = A^K^gap, see formula (lA.lip . and 



A„^(n) 



/3 





(61) 
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generating the d-metric g^/? = AtrA^gap (f58|) . 

If the metric and N-connection coefficients (I55p are stated to be those 
from an ansatz (p6]) (or ([M])), we should write °g = g®Nh (or °g = g®Nh) 
and say that the metric g = g ©iv h (150p was nonholonomicahy deformed in 
order to generate an exact solution of the Einstein equations for the canonical 
d-connection (or, in a restricted case, for the Levi Civita connection). In 
general, such metrics have very different geometrical and physical properties. 
Nevertheless, at least for some classes of 'small' nonsingular nonholonomic 
deformations, it is possible to preserve a similar physical interpretation by 
introducing small polarizations of metric coefficients and deformations of 
existing horizons, not changing the singular structure of curvature tensors. 
Explicit examples are constructed in Ref. [28j . 



4.2 The Geroch transforms as parametric nonholonomic de- 
formations 

We note that any metric °gap defining an exact solution of the vacuum 
Einstein equations can be represented in the form (j50p . Then, any met- 
ric °fi'o/3(^) (|B.2p from a family of new solutions generated by the first type 
parametric transform can be written as (j54p and related via certain polariza- 
tion functions of type (jSSp . in the parametric case depending on parameter 
0, i.e. rja{9) = {r]i{6),r]a{0)) and r]f{9)- Roughly speaking, any para- 
metric transform can be represented as a generalized class of anholonomic 
frame transforms additionally parametrized by and adapted to preserve 
the (n -|- m)-splitting structure|j The corresponding frame transforms (jB.SP 
and (IB.6j) are parametrized, respectively, by matrices of type (I60p and ()6ip . 
also "labelled" by 6. Such nonholonomic parametric deformations 



g 



5' 



'h :m = :m ©^(,) x^) 



are defined by the frame matrices. 







1i 

Oa~ 



(62) 



(63) 



/3 



generating the d-metric °ga/3 = °Acr °ggj3 and 



Ai^(n,( 



V\Va{u,9)\5a- 



(64) 



*It should be emphasized that such constructions are not trivial, for usual coordinate 
transforms, if at least one of the primary or target metrics is generic off-diagonal. 
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generating the d-metric °ga/3(^) = Aj^A^- °gaj3- Using the matrices ([63]) 
and (I64p . we can compute the matrix of parametric transforms 

B/=Ai^ ^A/, (65) 

Uke in (jB.Tp . but for "boldfaced' objects, where °A^' is inverse to °A^t, 
and define the target set of metrics in the form 

There are two substantial differences from the case of usual anholonomic 
frame transforms (j59p and the case of parametric deformations (j62p . The 
first one is that the metric g was not constrained to be an exact solution 
of the Einstein equations like it was required for °g. The second one is 
that even g can be restricted to be an exact vacuum solution, generated 
by a special type of deformations ([58]) . in order to get an ansatz of type 
(I44p . an arbitrary metric from a family of solutions °go/3(^) will not be 
parametrized in a form that the coefficients will satisfy the conditions (j45p . 
Nevertheless, even in such cases, we can consider additional nonholonomic 
frame transforms when g is transformed into an exact solution and any 
particular metric from the set | °Eapi0)} will be deformed into an exact 
solution defined by an ansatz (jUp with additional dependence on 9. 

By superpositions of nonholonomic deformations, we can parametrize a 
solution formally constructed following by the parametric method (from a 
primary solution depending on variables in the form 

°g(e) = e^(^''^''^) [e2 dx^ ® dx^ + eg dx^ ® dx^] (66) 
+64/10 [b* {x\v,9)]'^ 6v^6v 
+65 [b {x^ x^v, 9) - bo{x^,x^9)] ' 5/ ® 
Sv = dv + W2 {x'^ , x^ ,v,9) dx'^ + [x'^ jX"^ ,v,9) dx^ , 
6y^ = dy^ + n2 {x"^ , x^ , 9) dx'^ + ns [x'^ , x^ , 9) dx^ , 

with the coefficients restricted to solve the equations (fl5]l but depending 
additionally on parameter 9, 

e2r'{d) + e^i:"{9) = 0, (67) 
w'2{9)-wl{9) + w:,wl{9)-W2{9)wl{9) = 0, 

n'^{9)-nl{9) = 0, 

^"we use a "boldface" symbol because in this case the constructions are adapted to a 
(n + m)-splitting 



24 



for W2ie) = {b*{e))-^ {b{e)+bo{e))' and W3 = {b*{9))-^ {b{e)+bo{e)y, where, 



for instance, n'{9) = 92^3 (^) and 77-2 = d-in2{9). 

One should be noted that even, in general, any primary solution °g can 
not be parametrized as an ansatz (j44p . it is possible to define nonholonomic 
deformations to a such type generic off-diagonal ansatz °g or any g, defined 
by an ansatz (|50p . which in its turn can be transformed into a metric of type 
(|66p without dependence on 00 

Finally, we emphasize that in spite of the fact that both the paramet- 
ric and anholonomic frame transforms can be parametrized in very simi- 
lar forms by using frame transforms there is a criteria distinguishing one 
from another: For a "pure" parametric transform, the matrix BJ^ {u,9) 
and related Aq- and are generated by a solution of the Geroch equa- 

tions ([Rill . If the" pure" nonholonomic deformations, or their superposition 
with a parametric transform, are introduced into consideration, the matrix 
Acr(u) (IHT]) . or its generalization to a matrix Aa~ (f64l) . can be not derived 
only from solutions of ()B.4p . Such transforms define certain, in general, 
nonintegrable distributions related to new classes of Einstein equations. 

4.3 Two parameter transforms of nonholonomic solutions 

As a matter of principle, any first type parameter transform can be 
represented as a generalized anholonomic frame transform labelled by an 
additional parameter. It should be also noted that there are two possibilities 
to define superpositions of the parameter transforms and anholonomic frame 
deformations both resulting in new classes of exact solutions of the vacuum 
Einstein equations. In the first case, we start with a parameter transform 
and, in the second case, the anholonomic deformations are considered from 
the very beginning. The aim of this section is to examine such possibilities. 

Firstly, let us consider an exact vacuum solution °g (j44p in Einstein grav- 
ity generated following the anholonomic frame method. Even it is generic 
off-diagonal and depends on various types of integration functions and con- 
stants, it is obvious that it possess at least a Killing vector symmetry because 
the metric does not depend on variable . We can apply the first type pa- 
rameter transform to a such metric generated by anholonomic deforms (j59p . 
If we work in a coordinate base with the coefficients of °g defined in the 
form °g Q = °ggj3, we generate a set of exact solutions 



^^in our formuleis we shall not point dependencies on coordinate variables if that will 
not result in ambiguities 
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see ()B.2p . were the transforms ()B.7p . labelled by a parameter 9', are not 
adapted to a nonholonomic (n + m)-splitting. We can elaborate N~adapted 
constructions starting with an exact solution parametrized in the form (j54p . 

for instance, like °ga'i3' = ^^^^9gji , with AcT being of type ([6T]) with 
coefficients satisfying the conditions |58j). The target 'boldface' solutions 
are generated as transforms 



(68) 



where 



like in ()B.7p . but for "boldfaced' objects, the matrix °AJ^' is inverse to 

b' Ajb': a 



and there is considered the matrix 











^/\riaVa{0')\6a- 



where 'qi{u,6'),rja{u,6') and rjf{u,6') are gravitational polarizations of type 
(|55p l^ Here it should be emphasized that even °gai3{0') are exact solutions 
of the vacuum Einstein equations they can not be represented by ansatz 
of type (j66p . with 6 — >■ 6', because the mentioned polarizations were not 
constrained to be of type ([58]) and satisfy any conditions of type ([67|) F^ 

Now we prove that by using superpositions of nonholonomic and param- 
eter transforms we can generate two parameter families of solutions. This 
is possible, for instance, if the metric °ga'i3' form ([68]) . in its turn, was 
generated as an ansatz of type ([66]) . from another exact solution °ga"i3"- 
We write ^ ^ 

°ga'^'{e) = Bf{u,e) Bf{u,e) °g„";3" 

and define the superposition of transforms 

:gap{9',e) = B^-\e') Bf 



{9') B„? (9) B^^ (9) rg."/3'- 



(69) 



It can be considered an iteration procedure of nonholonomic parameter 
transforms of type (j69p when an exact vacuum solution of the Einstein 



^■^we do not summarize on repeating two indices if they both are of lower/ upper type 
^^As a matter of principle, we can deform nonholonomically any solution from the family 
^gap{d') to an ansatz of type (|66l) . 
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equations is related via a multi ^-parameters frame map with another pre- 
scribed vacuum solution. Using anholonomic deformations, one introduces 
(into chains of such transforms) certain classes of metrics which are not exact 
solutions but nonholonomically deformed from, or to, some exact solutions. 

4.4 Multi— parametric Einstein spaces and Ricci flows 

Let us denote by 6 = [^9 = 9', "^6, ...,9 = ^6,) a chain of nonholo- 
nomic parametric transforms (it can be more general as (j69p . beginning with 
an arbitrary metric ga"i3") resulting in a metric ga/si 9 ). Any step of non- 
holonomic parametric and/ or frame transforms are parametrize matrices of 
type ([Ml) . (f65]) or ([681) . Here, for simplicity, we consider two important ex- 
amples when go,^ ( 9 ) will generate solutions of the nonholonomic Einstein 
equations or Ricci flow equations. 



4.4.1 Example 1: 

We get a multi-parametric ansatz of type (llSp with /14 7^ and /15 ^ 
if ga/3( G ) is of type 

°g(V) = ei dx^ (g) dx^ e2fi'2(V,x*) dx^ (g) dx^ 
+£353 (V,x*) dx^ (g)dx^ 



ri 9 ,x\v 



+£5 



/ 



x\v]-fo{9,x') 5y'®5y\ 



,x^,v] 1 6v ® 6v 



(70) 



dv + Wk( 9 ,x\v] dx , 6y^ = dy^ + 



,x\v^ dx'', 



5v 



the indices with "hat" running the values = 2,3, where g^yu ,x 

are multi-parametric families of solutions of the 2D equation ([20p for given 
sources T4 ( 6* , x 



-hlCf,x') X 



? 1 9 ,x\vj = ?[o] [9 ,x 
T2(V,x^,?;)r I ,x\v 



x\v)-foi 9 ,x*) 



dv, 
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and the N-connection Nf = Wi{ 9 ,x^,v), Nf = ni{ 9 ,x^,v) computed 

Wi(9 , x\v) = }^ f , (71) 



Uk (V, x^ = nfc[i] (V, X*) + nk[2] (V, x^) x (72) 

(^^,x\v^ dv, 



,x\v 



f iy ,x\v) - Ux^) 



define an exact solution of the Einstein equations d?]). We emphasize that 
such solutions depend on an arbitrary nontrivial function 9 ,x'^,v] , 

with /* / 0, integration functions /o( 9 ,x*), /iq( 9 ,x*), <f[o] ^ 6* ,x* 

nfc[i] ^ ,x*^ , nfc[2] ^ and sources T2( 9 ,x^,w),T4 ^ 6* . Such 
values for the corresponding signatures = ±1 have to be defined by 
certain boundary conditions and physical considerations. We note that for- 
mulas ([7T]) and ([72|) state symbolically that at any intermediary step from 
the chain 9 one construct the solution following the respective formulas 
(I27p and (I28p . The final aim, is to get a set of metrics ()70p . parametrized 
by , when for fixed values of 0-parameters, we get solutions of type ([26]) . 
for the vacuum Einstein equations for the canonical d-connection. 

4.4.2 Example 2: 

< — > 

We consider a family of ansatz, labelled by a set of parameters 9 and x 

(as a matter of principle, we can identify the Ricci flow parameter x with any 
i — > 

9 from the set 9 considering that the evolution parameter is also related 
to the invariance of Killing equations, see Appendix [B]) . 

^g( V, x) = e^^"^'^''^''^) [ea dx^ ® dx^ + dx^ ® dx^] + (73) 
h4 9 ,x^ , x^ ,v,x^ 6v ® 5v 

Sv = dv + W2 9 jx"^ ,x^ ,vj dx^ + ^ 9 , x^, x^, u ^ dx^, 



= dy^ + 722 (|V,x^,x^,x) + ^2;^ 
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which for any fixed set 9 is of type (08]) with the coefficients are subjected 
to the conditions (j49p . in our case generahzed in the form 

e2V'"(V,x) + e3/(V,x) = A, 

hiCe)(t>Ce)/h^Ce)h^{7) = a, (74) 

«;^(V)-«;^(V) + u;3(VK(V)-«;2(VK(V) = 0, 

9 ,X) -n'iO ,x) = 0, 

for uh- = d^(l)/(p* , see ([25]) . being compatible with ([37|) and ([38]) and consid- 
ered that finally on solve the Einstein equations for given surces T4 = A 
and T2 = A. The metrics (|73p define self-consistent evolutions of a multi- 
parametric class of general solutions of the Ricci flow equations (|12p - (|14p 
transformed equivalently in the form (I20p - (l23p . The additional constraints 
(|74p define multi-parametric integral subvarieties (foliations) when the evo- 
lutions are selected for the Levi Civita connections. 



5 Summary and Discussion 

In this work, we have developed an unified geometric approach to con- 
structing exact solutions in gravity and Ricci flow theories following superpo- 
sitions of anholonomic frame deformations and multi-parametric transforms 
with Killing symmetries. 

The anholonomic frame method, proposed for generalized Finsler and 
Lagrange theories and restricted to the Einstein and string gravity, applies 
the formalism of nonholonomic frame deformations |20] [21] [23] [18] (see out- 
line of results in [19] and references therein) when the gravitational field 
equations transform into systems of nonlinear partial differential equations 
which can be integrated in general form. The new classes of solutions are 
defined by generic off-diagonal metrics depending on integration functions 
on one, two and three/ four variables (if we consider four or five dimensional, 
in brief, 5D or 4D, spacetimes). The important property of such solutions is 
that they can be generalized for effective cosmological constants induced by 
certain locally anisotropic matter field interactions, quantum fluctuations 
and/or string corrections and from Ricci flow theory. 

In general relativity, there is also a method elaborated in Refs. |31] [32] 
as a general scheme when one (two) parameter families of exact solutions 
are defined by any source-free solutions of Einstein's equations with one 
(two) Killing vector field(s) (for nonholonomic manifolds, we call such trans- 
forms to be one-, two- or multi-parameter nonholonomic deformations/ 
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transforms). A successive iteration procedure results in a class of solutions 
characterized by an infinite number of parameters for a non-Abelian group 
involving arbitrary functions on one variable. 

Both the parametric deformation techniques combined with nonholo- 
nomic transforms state a number of possibilities to construct "target" fami- 
lies of exact solutions and evolution scenarios starting with primary metrics 
not subjected to the conditions to solve the Einstein equations. The new 
classes of solutions depend on group like and flow parameters and on sets of 
integration functions and constants resulting from the procedure of integrat- 
ing systems of partial differential equations to which the field equations are 
reduced for certain off-diagonal metric ansatz and generalized connections. 
Constraining the integral varieties, for a corresponding subset of integration 
functions, the target solutions are determined to define Einstein spacetimes 
and their Ricci flow evolutions. In general, such configurations are nonholo- 
nomic but can constrained to define geometric evolutions for the Levi Civita 
connections. 

Acknowledgement: The work is performed during a visit at Fields 
Institute. 



A The Anholonomic Frame Method 



We outline the geometry of nonholonomic frame deformations and non- 
linear connection (N-connection) structures [181 US] • 

Let us consider a (n + m)-dimensional manifold V enabled with a pre- 
scribed frame structure ([1]) when frame transforms are linear on N^{u), 



e^{u) 

e\{u) Nl{n)e\{u) 
e\{u) 



(A.l) 
(A.2) 



where i, j, .. = 1, 2, n and a, b, ... = n + l,n + 2, ...n + m and u = {u" = 
{x^,y"')} are local coordinates. The geometric constructions will be adapted 
to a conventional n + m splitting stated by a set of coefficients N = {N^{u)} 
defining a nonlinear connection (N-connection) structure as a nonintergrable 
distribution 

TV =hYevV (A.3) 

with a conventional horizontal (h) subspace, KV, (with geometric objects 
labelled by "horizontal" indices and vertical (v) subspace v\' (with 
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geometric objects labelled by indices a,b, ...). The "boldfaced" symbols will 
be used to emphasize that certain spaces (geometric objects) are provided 
(adapted) with (to) a N-connection structure N. 

The transforms (jA.ip and ()A.2p define a N-adapted frame structure 

= (e^ = di- N^{u)da, ea = da) , (A.4) 

and the dual frame (coframe) structure 

e'^ = (e* = dx\ e" = dy" + Nl{u)dx') . (A.5) 

The frames (|A.5P satisfy the certain nonholonomy (equivalently, anholon- 
omy) relations of type ([2]), 

[e„, e/j] = BaBf^ - e^e„ = VFj^e^, (A. 6) 

with anholonomy coefficients 

Wt = daN^ and W^, = fi?, = e,(iVf ) - e,(Arf ). (A.7) 

A distribution ()A.3|) is integrable, i.e. V is a foliation, if and only if the co- 
efficients defined by N = {N^{u)} satisfy the condition 0°^- = 0. In general, 
a spacetime with prescribed nonholonomic splitting into h- and v-subspaces 
can be considered as a nonholonomic manifold [181 [T71 [S] • 
Let us consider a metric structure on V, 

g = g^^ (u) du^ ® du^ (A.8) 



defined by coefficients 

r n. . J- Aja Ajbh , Ajeh 1 

(A. 



Nfhe hah 



This metric is generic off-diagonal, i.e. it can not be diagonalized by any 
coordinate transforms if Nf{u) are any general functions. We can adapt 
the metric (jA.Sp to a N-connection structure N = {Nf{u)} induced by the 
off-diagonal coefficients in ()A.9P if we impose that the conditions 

g{ei, ea) = 0, equivalently, g.^ - N^hab = 0, 

where g^^ = g{d/dx^,d/dy"'), are satisfied for the corresponding local basis 
()A.4p . In this case A^^^^ = h"'^g. , where h""^ is inverse to hab, and we can write 
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the metric g ()A.9P in equivalent form, as a distinguished metric (d-metric) 
adapted to a N-connection structure, 

g = ga/3 {u) e" = g,j (n) ® + Kb (u) (g) e^ (A. 10) 

where gij = g(ej,ej) and hat = g(ea,eb) . The coefficients gap and g^^ = 
gaji are related by formulas 

ga/3 = A^A^-g^^, (A.ll) 

or 

5ii = ej-e/fl'y and /i^b = ej^e^-gab, 

where the frame transform is given by matrices (jA.ip with e-- = Sf and 
Co," = (^a~- We shall call some geometric objects, for instance, tensors, con- 
nections,..., to be distinguished by a N-connection structure, in brief, d- 
tensors, d-connections,... if they are stated by components computed with 
respect to N-adapted frames ()A.4p and (jA.Sp . In this case, the geometric 
constructions are elaborated in N-adapted form, i.e. they are adapted to 
the nonholonomic distribution (IA.3j) . 

Any vector field X = (hX, vX) on TV can be written in N-adapted 
form as a d~vector 

X =X"e„ = {hX = X'ei, vX = X^Ca). 

In a similar form, we can 'N-adapt' any tensor object and get a d-tensor. 

By definition, a d-connection is adapted to the distribution (IA.3j) and 
splits into h- and v-covariant derivatives, D = hD + vD, where hD = 
{Dfc = (t*^,L^^^} and vD = {Dc = (c'j^,C^^^} are correspondingly 
introduced as h- and v-parametrizations of the coefficients 

Lj, = (D,e,)Je\ L^, = (Dfceb)Je", C], = (D,e,)Je\ = (Dee,)Je«. 

The components V^^p ~ (-^jfe' -^fofc' ^jc ^6c) ' ^^^^ the coefficients defined 
with respect to ()A.5P and ()A.4p . completely define a d-connection D on a 
N-anholonomic manifold V. 

The simplest way to perform a local covariant calculus by applying d- 
connections is to use N-adapted differential forms and to introduce the d- 
connection 1-form F"^ = F'^^^e'*', when the N-adapted components of d- 
connection = (e^jD) are computed following formulas 

F^ .(^x) = (D„e^)Je^ (A.12) 
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where "J" denotes the interior product. We define in N-adapted form the 
torsion T = {T"} 

T° = De" = (ie° + r°^ Ae„, (A. 13) 

and curvature R = {7?."^} ([1]), 

7^"^ = TtTi = dTi- r^^ A r% . ( a. i4) 

The coefficients of torsion T ()A.13P of a d-connection D (in brief, d- 
torsion) are computed with respect to N-adapted frames (jA.Sp and (|A.4p . 

rpi ri ri rpi rpl /^i rpa r\a 

jk ^ jk kji ^ ja ^ aj ^ ja) ^ ji ji) 

dN"- 

bi — ib — "oyb ^ bi^ ^ be — bc~ '-^ cb^ [A. 10) 

where, for instance, T^^ and T'^^ are respectively the coefficients of the 
/i(/i/i)~torsion hT{hX,hY) and v (wj-torsion vT{ vX, vY). In a similar 
form, we can compute the coefficients of a curvature R, d-curvatures. 

There is a preferred, canonical d-connection structure, D, on a N- 
anholonomic manifold V constructed only from the metric and N-connection 
coefficients [gij,hab,N^] and satisfying the conditions Dg = and T^jf^ = 
and T^^ = 0. It should be noted that, in general, the components T^ja; 
and T'^- are not zero. This is an anholonomic frame (equivalently, off- 
diagonal metric) effect. Hereafter, we consider only geometric constructions 
with the canonical d-connection which allow, for simplicity, to omit "hats" 
on d-objects. We can verify by straightforward calculations that the linear 
connection F'^^^ = (^L^jk^ L^^k^ ^jc^ ^bcj ^^^^ the coefficients defined 

Defe(ej) = LjfcBi, De,(efe) = L^^e^, De,(ej) = Cj^e,, De^(efc) = C^^Ca, 
where 

= 1^9" i^kgjr + ejgkr - ^rQjk) , 

L^fc = eb{N^) + hi'''{ekH,-hdcebNt-hdbecNt), (A.16) 

1 1 
■' Z z 

uniquely solve the conditions stated for the canonical d-connection. 
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The Levi Civita hnear connection V = { i^^^}, uniquely defined by the 
conditions = and \/g = 0, is not adapted to the distribution (|A.3p . 
Denoting ,r°^ = ( ,L}„ ,L^, ,q„ for 

Vefc(ej) = |L*^ej + |L";,ea, Vefe(ef,) = i-^^L^i + i-^bfe^a, 

Vei,(ej) = iCjf^ej + |C";,ea, Vec(efe) = ,Cl^ei + .C^cea, 
after a straightforward calculus we get 

.Ltk = Ll^ + \{6t5'd + hcdh'''')[Llk-eb{m)]. 

= -3-l{[Ll^-ea{N^)\K,+ [Ll^-e,{N-)\h,a], 

where r2|^ are computed as in the second formula in ()A.7p . 

For our purposes, it is important to state the conditions when both 
the Levi Civita connection and the canonical d-connection may be defined 
by the same set of coefficients with respect to a fixed frame of reference. 
Following formulas ()A.16P and ()A.17p . we obtain equality iF^^ = V^^p if 

^% = (A.18) 

(there are satisfied the integrability conditions and our manifold admits a 
foliation structure), 

.Clb = CL = (A.19) 

and LJjj — ea{Nj) = 0, which, following the second formula in ()A.16p . is 
equivalent to 

ekhc - hdc ebNi - h^b e^N^ = 0. (A.20) 

We conclude this section with the remark that if the conditions (jA.lSp . 
([Al9]) and (1X20]) hold true for the metric (|X8]) . equivalently ([AlO]) . the 
torsion coefficients ()A.15P vanish. This results in respective equalities of the 
coefficients of the Riemann, Ricci and Einstein tensors (the conditions ([6]) 
being satisfied) for two different linear connections. 
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B The Killing Vectors Formalism 



The first parametric method (on holonomic (pseudo) Riemannian spaces, 
it is also called the Geroch method [31]) proposes a scheme of constructing 
a one-parameter family of vacuum exact solutions (labelled by tilde and 
depending on a real parameter 9) 

':m= ^gape'^^e^ (B.l) 

beginning with any source-free solution °g = { "Qa/s} with Killing vector 
^ = {^a} symmetry satisfying the conditions ,£ = (Einstein equations) 
and V^( °g) = (Killing equations). We denote this 'primary' spacetime 
"g^ia) and follow the conventions: The class of metrics °g is generated 
by the transforms 

°5a/3 = 4°' in, e) Bf {U, 9) °(7a'/3' (B.2) 
where the matrix is parametrized in the form when 

°gaP = AA^^( °ga(3 - X'^iaip) + (B.3) 

for 

A = \[{cos9 -ujsiii9f + \^s\n^9]-^ 
fir = A~^^r + cut sin 29 — j3r sin^ 9. 

A rigorous proof [3T| states that the metrics ()B.ip define also exact vacuum 
solutions with = if and only if the values ^^a, ar, /^r from (|B.3p . subjected 
to the conditions A = ia^p OJ = a-y , ^'^ fx^ = A^ + o;^ — 1, solve the 

equations 

Vau = e^p^rf V^r, V[„a^] = ^e^p^r V^r, (B.4) 

where the Levi Civita connection V is defined by °g and ea/s^yr is the abso- 
lutely antisymmetric tensor. The existence of solutions for ()B.4p (Geroch's 
equations) is guaranteed by the Einstein's and Killing equations. 

The first type of parametric transforms ()B.2p can be parametrized by 
a matrix B^ with the coefficients depending functionally on solutions for 
()B.4p . Fixing a signature g^^ = diag[zizl, ±1, .... it 1] and a local coordinate 
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system on {V, °g,S,a), we define a local frame of reference e^' = A^{u)da, 
like in ([T|), for which 

°ga'f3' = A^A^gap. (B.5) 

We note that A^ have to be constructed as a solution of a system of 
quadratic algebraic equations (|B.5P for given values gaj3 and °ga'p'- In a 
similar form, we can write Cq, = Acr(^, "0)190, and 

The method guarantees that the family of spacetimes {V, °g) is also vac- 
uum Einstein but for the corresponding families of Levi Civita connections 
V. In explicit form, the matrix B^^ {u,9) of parametric transforms can be 
computed by introducing the relations (jB.Sp . (|B.6P into ()B.2p . 

B^^' = A^ A^ (B.7) 

where A^' is inverse to A^. 

The second parametric method |32j was similarly developed which yields 
a family of new exact solutions involving two arbitrary functions on one vari- 
ables, beginning with any two commuting Killing fields for which a certain 
pair of constants vanish (for instance, the exterior field of a rotating star). 
By successive iterating such parametric transforms, one generates a class 
of exact solutions characterized by an infinite number of parameters and 
involving arbitrary functions. For simplicity, in this work we shall consider 
only a nonholonomic version of the first parametric method. 
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